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We study dynamical properties of systems with many interacting Fermi-particles under the 
influence of static imperfections. Main attention is payed to the time dependence of the Shannon 
entropy of wave packets, and to the fidelity of the dynamics. Our question is how the entropy and 
fidelity are sensitive to the noise. In our study, we use both random matrix models with two-body 
interaction and dynamical models of a quantum computation. Numerical data are compared with 
analytical predictions. 

O ' I. INTRODUCTION 

o : 

(N ■ 

^ , There are two sources for the appearance of statistical properties in quantum isolated systems. The first one which 

■ serves as a base for the traditional statistical mechanics, is the thermodynamic limit in which the number of particles is 
, infinitely large. In this case even integrable systems may be treated statistically. The origin of this effect is an infinite 

■ number of independent frequencies in the dynamics of a system. It is clear that in this case any small interaction 
, between particles, or with a heat bath, gives rise to strong statistical properties. 

■ Another situation occurs when the number of particles is finite and not very large. In this case the role of interaction 
is crucial, and chaos arises under some conditions. For systems with the well defined classical limit, this chaos is known 
as quantum chaos and can be compared with the classical one. This situation is known in literature as many-body 
chaos, in contrast to the thoroughly studied case of one-body chaos which refers to a single particle in an external 

, field. The well known example is the hydrogen atom in a strong magnetic field. 
' In many-body systems the density of many-particle energy levels increases extremely fast (typically, exponentially) , 
CO both with an increase of the number of particles and excitation energy. For this reason, the interaction between 
' particles can lead to a strong mixing between many-particle basis states, thus resulting in chaotic eigenstates. The 
"nI ' latter term refers to the fact that the components of such eigenstates can be practically treated as pseudo-random 
variables, this leads to strong statistical properties such as the relaxation of the system to a steady-state distribution. 
f| Typical examples of such systems are compound nuclei, complex atoms, atomic clusters, isolated quantum dots, etc. 
Qh. Recent calculations for complex atoms [1], multicharged ions [2], nuclei [3], Bose-Einstein condensates [4], and spin 
4!) ■ systems [5,6] have confirmed the dynamical origin of statistical laws in isolated systems (see details in [7,8] and the 
^ [ review [9]). 

2 ' The onset of chaos for highly excited states and for many-particle spectra has been thoroughly studied by making 
qh] use of the Two-Body Random Interaction (TBRI) model which was invented long ago [10]. In this model all two-body 
• • . matrix elements are assumed to be independent and random variables, however, the model is essentially different 
from standard random matrix models where the two-body nature of interaction is not taken into account. One of the 
important results obtained in the frame of this model (see, for example, [11] and references therein), is the Andcrson- 
^ ■ like transition which occurs in the space determined by many-particle states of the unperturbed part Hq. Due to a 
5^ ' random character of the perturbation, this delocalization transition can be treated as the transition to a global chaos. 
Above some critical value of the perturbation, the number Npc of principal components in eigenstates is typically 
large. For example, Npc ~ 150 in the Ce atom [1] and Npc w 10'' — 10^ in heavy nuclei [3]. 

Recently, the theory of many-body chaos has been extended to quantum computers. One can naturally expect 
[12,13] that due to a very high density of energy levels, any kind of perturbation may lead to decoherence effects 
thus destroying the quantum computation. For this reason it is very important to search for the conditions when 
the role of chaos can be significantly reduced [14,15]. So far, the study of the many-body chaos has been mainly 
restricted by statistical properties of the energy spectra and eigenstates. On the other hand, in view of experimental 
applications, one needs to know what are the dynamical properties of quantum systems with strongly interacting 
particles. Especially, this is very important for quantum computation [16], for which quantum protocols assume a 
high stability of a long-time dynamics under the influence of an environment or any kind of static imperfections. 
Below, we analyze how the static disorder influences the dynamical properties of quantum systems with complex 
behavior. 
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II. DYNAMICS OF SYSTEMS WITH TWO-BODY INTERACTION 



In what follows we assume that the models under consideration can be represented by the Hamiltonian separated 
in two parts, 

H = Ho + V (1) 
where Hq describes the "unperturbed" part and V stands for the "perturbation" . Then V can be represented as 

Vik = {l\V\k), (2) 

where |/) and \k) stand for basis states of Hq (or, "unperturbed" states). Correspondingly, exact states \a) of the total 
Hamiltonian H are expressed as 

\a)=C'^\k). (3) 

The coefficients give the expansion of an exact state in terms of the basis states (for a fixed), or the expansion of 
a basis state in terms of the exact states (for k fixed). 

In the case when the Hamiltonian H corresponds to a system of interacting Fermi-particlcs, basis states can be 
constructed as |fc) — a^^...a^ |0) from the ground state |0) with a'^ as the creation operator. In application to 
quantum computer models the Hamiltonian Hq typically describes a number of non-interacting qubits, and V stands 
for the inter-qubit interaction needed for quantum computation (for other variants, sec below). In this case the basis 
state |fc) is a product of single qubit states, with af as the spin-raising operator (if the ground state |0) corresponds 
to all spins down). 

In principle, the knowledge of the state matrix and of the corresponding energy spectrum E"" gives a complete 
information about the system. In particular, if an initial state ^(0) is some basis state \ko), the evolution of the 
'J'— function is described by the expression, 

m = E ^nC^ eM-iE'^t). (4) 

Here and below we assume that fi = l. As one can see, the probability 

Wm = \Amf = I (m|*(i)) I' (5) 
to find the system at time t in the state |m) is determined by the amplitude 

= {m\eM-^Ht)H = ^ [C^l^ exp(-zi^"t). (6) 

a 

When the perturbation is weak, only the component Am with m = fco is large, therefore, the evolution of the system 

is close to the periodic one. However, with an increase of the perturbation strength, the number of components Am 
with relatively large amplitudes increases, and the dynamics of wave packets can be quite complicated. In order to 
characterize the dynamics of a system in the unperturbed basis |A:), one uses different quantities. The main interest 
is to known how many unperturbed states are involved in the dynamics depending on time. One of the commonly 
used quantities is the Shannon entropy of the packet, 

S{t) = -^Wm^I^Wm = -WolnWo - ^ Wm^riWm; ^ = 1. (7) 
m mj^O rn 

Here we introduced the return probability Wo(t) = |^o(i)|^ for the system to be in the initial state |A:o), due its special 
role in global properties of the dynamics. One can see that the effective number Np of unperturbed states involved in 

the dynamics can be estimated as Np w In 5(f). 

Another quantity which nowadays is under the close investigation, is the so-called fidelity, 

m = \{%imu{m'' (8) 

where 4'„(f) stands for the '5— function of the unperturbed system, and ^p{t) for that of the perturbed one. In 
the case when perturbation is weak, this quantity serves as a measure of the sensitivity of the system to a given 
perturbation [17]. There are many results obtained for the fidelity, both analytical and numerical ones, obtained for 
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different systems (sec, for example, [18] and references therein). The main interest is whether the fidelity can serve 
as a good indicator of the quantum chaos. Note, that the fidelity depends both on the type of the perturbation, and 
on the form of an initial packet. In what follows, we discuss this quantity in application to systems of interacting 
Fermi-particles and to a model of quantum computation. 

By comparing the definition (8) with Eqs.(5) and (6), one can see that in the case when the initial state is a specific 
basis state |fco), the fidelity is nothing but the return probability [19], 

T{t) = VFo(i) = 1^0 Wl' ; Ao{t) = J2 |Cfc"J'exp(-zi^"t) « / Pk,{E) exp{-iEt)dE. (9) 

Here we replaced the summation by the integration that can be done when the number of components C"^ is large 

for a fixed fco. Generically, these components fluctuate strongly, therefore, they can be treated as pseudo-random 
quantities. In fact, this condition of a large number of pseudo-random components in exact eigenstatcs can be used 
as the definition of chaos in quantum systems. In this case, the time dependence of Wo(t) is entirely determined by 
the Fourier representation of Pkg{E) = P{E,Ekg) where E is the energy of exact eigenstates and Ek^ is the energy 
corresponding to the unperturbed state |A;o) . This quantity is known in the literature as the strength function (SF), 
or, as the local spectral density of states, 



PiE,Ek,)^\Cl^>\^p{E). (10) 

Here p{E) is the density of states of the total Hamiltonian H, and the average is performed over a number of states 
with energies close to E. In next Section we discuss main properties of the fidelity (9) in application to close system 
of interacting Fermi-particles [19]. 

III. INTERACTING FERMI-PARTICLES 

A. Two-body random interaction model 

Let us assume that our model consists of finite number Np of Fermi-particles that occupy M single-particle states 
\s) characterized by the corresponding energies eg. Then the unperturbed Hamiltonian Hq describes non-interacting 

particles, and the interaction between particles is represented by V. In this case the basis many-body states can 
be constructed by the Slater determinant, \k) = al^ . . . al^ |0), and the Hamiltonian takes the form 

M 

Ho = ^e^ajas. (11) 

Correspondingly, the interaction term can be represented in this basis as 

V =^^Vs,s2S3s^aiyg^as^as^ (12) 

where al^ and ag^. are the creation-annihilation operators. 

The total Hamiltonian H = Hq + V with (11) and (12) appears in many physical applications such as complex 
atoms, nuclei, atomic clusters etc. In fact, the form of H discussed above is known as the mean field approximation 
for complex quantum systems of interacting particles. In this description, the unperturbed part Hq represents the 
zero-order mean field for the excited states with the ground state Ex , and the residual two-body interaction is given 
by V . Therefore, the single-particle energies in such applications are, in fact, renormalized quasi-particle energies 
(sec details, for example, in Ref. [21]). The considered model turns out to be very useful for understanding of the 
onset of many-body chaos due to a two-body interaction. 

For the analytical analysis one assumes that all two-body matrix elements Vg^s^s^si are random Gaussian numbers 
with the zero mean and the variance V^^ = (^51528354)- Although this assumption of the randomness is quite strong, 
however, it was found that many of the properties of this TBRI model are quite generic and occur in different 
dynamical (without any random parameters) systems. The results obtained in the frame of this model can be also 
extended to interacting Bose-particles, see, for example, Ref. [20] 

Without the loss of generality one can assume that the single-particle spectrum is non-degenerate, with the constant 
mean level spacing do = {^s+i — ^s) = 1- Therefore, the model is defined by four parameters, M, Np, do and Vq. 



3 



Note that the number N of many-body states increases very strongly with an increase of number of particles Np and 
number M of single-particles states, N = jv [(m-n )! • ^'^^ ^^^^ reason even for a relatively small number of particles 
the size of the Hamiltonian matrix is large, and exact eigenstates may consist of many unperturbed basis states, thus 
providing us with a possibility to describe the system statistically. 

Due to a two-body nature of the interaction, the Hamiltonian matrix has features that are different from those 
of the standard random matrices that are typically used when studying chaos in quantum systems. Specifically, the 
_ff— matrix is sparse, the off-diagonal matrix elements of V in the representation of many-particle basis states are not 
completely independent, and the global structure is a band- like (see details in Refs. [22,8]). This makes the analysis 
more complicated, on the other side, such matrices give a better description of real physical systems, in comparison 
with standard random matrix ensembles. 

One of the important results obtained recently in the frame of this model is the Anderson-like transition which 
occurs in the Fock space determined by many-particle states of Hq (see, for example, [11] and references therein). 
The critical value Vcr for this transition is determined by the density of states pf = dj^ oi those basis states which 
are directly coupled by the two-body interaction. When the interaction is very weak, Vq <C rf/ , exact eigenstates are 
very close to the unperturbed ones, consisting of a small number of basis states. With an increase of the interaction 
the number Npc of principal basis components in exact eigenstates increases, and for Vq > df the eigenstates may be 
treated as pseudo-random ones. In this case the number Npc can be estimated as Npc '^T/D where F is the spreading 
width of the strength fiuiction (10). In such chaotic eigenstates any external weak perturbation is exponentially 
enhanced. The enhancement factor can be estimated as ^Npc oc see Ref. [23] This huge enhancement have 

been observed in experiments when studying the parity violation in compound nuclei (see, for example, review [24]). 
In what follows, we discuss the dynamics of the model above this threshold of chaos. 



B. Fidelity 

We start now with the expression for the fidelity (9) at small times. According to the perturbation theory, one can 
easily get, 

Wo(i) « 1 - Ali^ (13) 
where is the width of the strength function (SF) in the energy space determined as 

^l=Y.^lk,- (14) 

The second moment A|; of the SF for the TBRI model can be found explicitly [7] , 

A| = \v^Np{Np - 1)(M - Np){M - JVp + 3) (15) 

where is the variance of the off-diagonal matrix elements of the two-body interaction V. It is interesting that for 

Fermi-particles the width A^; turns out to be independent of a specific basis state jfco). 

The expressions (13) and (14) are universal, they are correct for any kind of the perturbation V . Thus, the quadratic 
time dependence of the fidelity on a small time scale occurs independently on whether the dynamics is chaotic or 
regular. This is understandable since the difference between periodic and chaotic motion can be detected on a large 
time scale only (in principle, for t — s- oo). 

The situation is much more complicated for large times. In order to understand how strong the decrease of the 
fidelity on a large time scale, one needs to know the form of the strength function in the energy space. For a long 
time it was assumed that the SF has generically the Breit-Wigner form (or, the same, the Lorentzian form). On the 
other hand, recently it was understood that in many physical applications the form of the SF can be quite close to 
the Gaussian form (see, e.g. [25]). This is because in contrast to full random matrices, physical interactions V always 
have a finite width in the basis of Hq. This point was the reason for Wigner to introduce the so-called Band Random 
Matrices (BRM) with a finite width. Therefore, one more control parameter arises which is the energy associated 
with the width of perturbation, in addition to the width of the Lorentzian (see, for example, Ref. [25] and references 
therein). Full analytical treatment of the form of the SF for the BRM is given in Ref. [26] Specifically, it was shown 
that with an increase of the perturbation the form of the SF changes from the Lorentzian to the Gaussian (apart from 
extremely strong perturbation, when the non-physical semicircle form arises). 
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Similar transition from the Lorentzian to the Gaussian has been analytically found for the TBRI model discussed 
above [27], although the analytical expression in the closed form is not known. In order to evaluate the fidelity, in 
Ref. [19] the following phenomenological expression has been suggested, 



P{E) = B 



cxp 



{E - Eof + r2/4 



(16) 



where Eq is the energy that corresponds to the basis state jfco). Strictly speaking, this formula is valid near the center 
of the energy spectrum that has the Gaussian form, otherwise one should take into account additional distortion 
effects [7]. Due to the normalization conditions, / P{E)dE = 1 and / E'^P(E)dE = A|;, only one of three parameters 
B, r and a is free. The relations between these parameters are: 



1 

B 



1 _ $ 



r 



r V 8a2 



and 



(17) 



(18) 



where ^{z) is the error function. 

In the case of a relatively small (but non-perturbative) interaction, the form of P{E) is close to the Lorentzian. In 

this case F <C cr and T plays the role of the half- width Fq of the Lorentzian, F « Fq. On the other hand, for a strong 
perturbation the SF has the Gaussian form with a Ae and formally F ^ u. The value of Fq for the TBRI model 
is determined by the Fermi golden rule, 



ro{E)^2n\Vk,j\'pf{E) 



(19) 



where Pf{E) is the density of states directly coupled to the basis state jfco) by the two-body interaction V. 

In the case of a relatively weak interaction (the Lorentzian form of the SF) for large time, t tc = l/c, the 
evaluation of the fidelity gives, 



1 r2 


2 
E 



Wo{t) « exp ( - Fot ) « exp {-Tot) 



(20) 



where the correction term - 



—^T= is small. 



In the other limit case of a strong interaction (the Gaussian form of the SF), F a w A^;, we have tc ^ ^■ 
Therefore, the leading dependence of Wo{t) is the Gaussian, 



Wo{t) ~ exp(-A^i^), 
and only for a very long time ^ it becomes the exponential function [19] , 



(21) 



Woit) 



^2^2 /J r2 



SAl^'HiAl 



Tt 



(22) 



It is important to note that in this case the return probability Wo{t) has large correction factor exp (l^), 
addition to the standard decay law Wo{t) = exp(— Ff). 



C. Shannon entropy 



For the Shannon entropy (7) one needs to know the time dependence of all amplitudes Wm{t), see Eq.(5). In order 
to find Wm{t), in Ref. [25] the theory has been developed which is based on the representation of the dynamics of the 
system as a flow in the Fock space of many-body states. This approach can be compared with that developed in Ref. 
[28] where the dynamics in the many-particle basis was represented by a flow on the Caley tree. For the exponential 
dependence Wo{t) = exp(— Ft), the solution for W„(t) has simple form. 
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Wn{t) = ^W-o(t). (23) 
n! 

As one can sec, the dynamics of the system can be expressed in terms of the return probabiUty VFo(t) whieh is a 
particular case of the fidehty. Note, that the quantity F stands here for the effective width of the SF, and can be 
either the half-width of the Lorentzian (F = Fq), or the square root of the variance of the SF, (F = A^), depending 
on its form (Lorentzian or Gaussian). 

As a result, the expression for the Shannon entropy reads, 

S{t)f^n\nNf+rt-e-^'y^\n^, (24) 

where Nf is the number of basis states directly coupled to the basis state |A;o). Two last terms in the right-hand-side 
of Eq.(24 ) turn out to be smaller than the first one, therefore, one can write, 

S{t)^TtlnNf{l + f{t)}!vTtlnNf (25) 

with some function /(t) <C 1 which slowly depends on time. In this estimate for the increase of entropy, the influence 
of fluctuations of Wm is not taken into account. It can be shown that for the gaussian fluctuations of the coefficients 
Am with the variance given by their mean-square values, and Npc{t) « exp(5(i)) » 1, the entropy has to be corrected 
by a small factor of the order of In 2 (sec, for example, [29]). 

One should note that at small times the expression for the entropy has to be modified since on this time scale the 
function Wo{t) has always the form Wo{t) = exp(— A|, t^). To do this, one needs to replace Tt in Eq.(23) by a more 
accurate expression, — In(Wo). This gives A'^t'^ for small times, t <§; F/A^, and Tt for large times, t ^ r/A|;. 
Therefore, at small times, t <C F/A|, the entropy is given by the expression [25], 

Sit)^ Alt^ {l + 0{\n{t^))). (26) 

As one can see, the "exact" expression (24) for the entropy is quite complicated. Instead, in Ref. [25] it was proposed 
to use the following simple expression which gives good approximation for systems with a small number of particles, 

S{t) = -Wo{t)lnWo{t) - (1 - Wo(i))ln ( ^^~JZ ) (27) 

\ Npc J 

Here NpT'^ is the maximal value of the number of principal components after the saturation of the entropy to its 
maximal value. Quite often, this phenomenological parameter can be estimated independently. 



IV. THE MODEL OF QUANTUM COMPUTATION 
A. Description of the model 

Let us now apply the above analysis of the fidelity and entropy to the model of quantum computation. This model 
describes a onc-dimensional chain of L identical 1/2-spins placed in an external magnetic field. For a selective resonant 
excitation of spins, it is suggested to choose the time independent part — B^{x) of a magnetic field with the constant 
gradient along the a;-direction. This provides different Larmor frequencies for different spins, Wfe = jnB^ = too + ak, 
where 7„ is the spin gyromagnetic ratio and ak = Xk- One can arrange relative directions of the chain with respect to 
the 2;-axis in a way that the dipole-dipole interaction is suppressed, and the main interaction between nuclear spins 
is due to the Ising-like interaction mediated by the chemical bonds. This model was first proposed in [30] (see also 
[31-33]) as a simple realization of a solid-state quantum computation. 

In order to implement a time-dependent quantum protocol, the spins are assumed to be subjected to a transversal 
circular polarized magnetic field. Therefore, the total magnetic field has the following form [34,32,33], 

B{t) = [6^ cos(V + ipp), -6^ sm{upt + ipp), B'{x)]. (28) 

In the above expression, If^, Vp, and (fp are the amplitudes, frequencies and phases of a circular polarized magnetic 

field, respectively. The latter is given by a sum of p = 1, P rectangular time-dependent pulses of length fp_|-i — tp, 
rotating in the {x, y)— plane, in order to make selective excitations of specific spins. Thus, the quantum Hamiltonian 
of the system has the form. 
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L-1 



fe=0 L n>k 



p=l fe=0 



where the "pulse function" 0p(t) equals 1 during the p-th pulse, for tp < t < tp+i, otherwise, it is zero. The 
quantities Jk,n stand for the Ising interaction between two qubits , cok are the frequencies of spin precession in the 
B^— magnetic field, ftp is the Rabi frequency of the p-th pulse, /^'^'^ = (l/2)o-^'^'^ with cr^'"'^ as the Pauli matrices, 
and Ik=Ik± ill 

For a specific j3-th pulse, it is convenient to represent the Hamiltonian (29) in the coordinate system that rotates 
with the frequency Up. Therefore, for the time tp < t < tp+i of the p-th pulse our model can be reduced to the 
stationary Hamiltonian, 

n(p) = - EiCkll + all K) - 2 E (30) 

fc=0 n>k 

where = (w^ — Vp), a = fipcos^p, and /3 = fipsincpp. 

Below we restrict our considerations by the Hamiltonian for a single pulse, by choosing tpp = 0. We also assume a 
constant interaction between nearest qubits, Jk,n = J5k,k+i- Then the Hamiltonian takes the form, 

L-1 L-1 

^^"^ = E [ - ^'^^fc - '^Jikii+i] - E ^p^k- (31) 

fe=0 fe=0 

In the ^-representation the Hamiltonian matrix of size 2^ is diagonal for Op = 0. Therefore, below we define Hq 
and V as 

L-1 L-1 

no = Y.[-^kil-2Jiiii^,] V = -J2^Ik (32) 

fe=o fe=o 

where we omitted the index p. The unperturbed basis (in which Tip is diagonal) is reordered according to an increase 
of the index s which is written in the binary representation, s = iL-i,iL-2, ■■■,io (with = or 1, depending on 
whether the single-particle state of the i— th qubit is the ground state or the excited one). The parameter fl thus is 
responsible for a non-diagonal coupling, determining matrix elements of V as Vkn = Vnk = — ^0/2 with n ^ k. As one 
can see, in contrast with the TBRI model discussed above, in the z— representation the interaction between particles 
is absorbed by Ho, and V describes the coupling to the external magnetic field. 

There are two regions of the parameters of physical interest. The first one is known as the so-called non-selective 
regime which is defined by the conditions, O ^ 5ujk ^ J ■ This inequality provides the simplest way to prepare a 
homogeneous superposition of 2^ states needed for the implementation of both Shor and Grover algorithms. The 
analytical and numerical treatment of the model (32) in this regime has shown [34] that constant gradient magnetic 
field (with the non-zero value of a) strongly reduces the effects of quantum chaos. Namely, the chaos border turns out 
to be independent of the number L of qubits, in contrast to the models thoroughly studied in Ref. [13] In particular, it 
was shown that quantum chaos occurs for a very large coupling O = 100 and very strong interaction J = 100 between 
qubits. For these parameters the Wigner-Dyson distribution between nearest energy levels has been observed, the fact 
that can serve as a numerical proof of the quantum chaos. The transition to chaos for = 100 in dependence on J 
can be analytically understood with the use of the transformation to the "mean field basis" in which the Hamiltonian 
is diagonal in the absence of the inter-qubit interaction, J = 0. In this basis the term with J ^ plays the role of 
the interaction between L blocks that correspond to quantum numbers, and the structure of the Hamiltonian matrix 
is similar to that of the TBRI model, sec details in Ref [34]. 

Another region characterized by inequalities ftp <C Jfc_„ <C a <C Wfc, is referred as the selective excitation, see Ref. 
[32] In this regime each magnetic pulse acts selectively on a chosen qubit, resulting in a resonant transition. According 
to the quantum protocol, many such resonant transitions take place for different p pulses, with different values of 
Vp ^ LUk- The detailed analytical analysis [35] has revealed that in this regime the perturbation theory works very well 
for many pulses, thus, indicating that there is no any effect of the quantum chaos. Therefore, the implementation of 
the constant gradient magnetic field is very effective in reducing any kind of decoherence. 

B. Shannon entropy and fidelity 

Now we discuss numerical data for the above model of quantum computation, by paying the main attention to the 
time-dependence S{t) of the Shannon entropy and fidelity J^{t). As we have discussed, the time-dependence of the 



7 



Shannon entropy provides us with the information about the cvohition of the wave packets in the unperturbed basis. 
Since the latter is determined by the Hamiltonian TYo in Eq.(32), the time-dependence of the entropy is directly related 
to the O— dependent term. The main interest is in the role of the interaction J, as well as of static random terms which 
we add to Specifically, each of non-zero off-diagonal elements is now the sum of two terms, => ilo+^p where 
stands for random variables with the zero mean and variance = (^p). This kind of randomness corresponds to the 
imperfections in the Rabi frequencies. 




FIG. 1. Normalized Shannon entropy for L = 8, o = l,f2o = 100, and different values of J: (a) J = 0, (b) J = 10, (c) 
J = 100. Curves with squares correspond to the perturbation ap = 5, with circles to ap = 10, and with triangles to <Tp = 20. 

In the first line, we discuss how the Shannon entropy depends on time if the initial state is the basis state chosen, 
for simplicity, at the center of the energy spectrum. Typical dependencies for S{t) are shown in Fig.l where we have 
normalized the entropy to its maximal value Smax determined by the total number of many-body states. The value 
Smax corresponds to the case when all basis states are equally populated, with the standard normalization that the 
total probability is one. As one can see from the data, the entropy saturates to its maximal value in two cases. The 
first one occurs for the dynamical chaos which is due to a large value of the inter-cubit interaction, J = 100. In this 
case there is no difference whether we have disorder ap or not. Another case is when the disorder is strong, ap = 20. 
In this situation, the saturation occurs even for J = 0. These results are quite instructive since they demonstrate 
equivalence of the dynamical chaos to a disorder. It should be pointed out that the entropy does not reach the 
maximal value S{t) = 1 (in the saturation) which may be explained by a suppressed chaos for the states with the 
energies close to the edges of the spectrum. The most important result is that for strong chaos (dynamical or due to 
the disorder), the system can be well described by statistical methods. 

The dept Smin of the first minimum in the time dependence of S{t) can be used as a relative measure of chaos in 
the system. In Fig. 2 the ratio Smin/ Smax is plotted against J, thus showing the transition to chaos. Surprisingly, 
the transition turns out to be quite smooth, in contrast to the transition measured in terms of an effective number of 
components in exact eigenstates [34]. The latter transition has revealed two borders, one is due to the "delocalization" 
of the eigenstates in the imperturbed basis (for J > 15), and the second one which is due to the quantum chaos in 
the eigenstates (for J ~ 100, when the level spacing distribution has the Wigner-Dyson form), see details in Ref [34]. 
Therefore, the Shannon entropy can serve as the indicator of the delocalization, rather than of the quantum chaos. 
Indeed, one can have a good relaxation even in a completely integrable system, if one uses the basis which is "very 
far" from the basis corresponding to the total Hamiltonian. 

Let us now analyze how good is the correspondence of the above data to the analytical expression (25) obtained 
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for the TBRI model. This expression gives the Hnear increase of the entropy before the saturation. The value of F in 
Eq.(25) can be expressed via the variance (14), since in our case V is formally large = 100). Taking into accoimt 
that in each row of the matrix V there are only Nf = L non-zero elements, one can easily compare the expression (25) 
with the numerical data, see Fig. 3. One can see that when the disorder is absent or small, ap = 0; 15, the linear slope 
of S{t) is slightly different from that predicted by Eq.(25). On the other hand, when the disorder is strong, ap = 50, 
the correspondence between the analytical expression and numerical data is quite good. Thus, the theory developed 
in the frame of the TBRI model appears to be also valid for the models that are not fully random. Similar effect was 
also observed in the Bose-Einstein condensate model [20] . Namely, for a relatively strong interaction between bosons, 
the Shannon entropy of the wave packet which was initially in the condensate, was found to increase linearly with 
time in a good agreement with the analytical expression (25). This proves an effectiveness of the approach developed 
in Ref. [34] in application to both dynamical and random systems. 

1 , , , , , , 




FIG. 3. Small time scale for the Shannon entropy where the linear increase is expected before the saturation. The data 
are given for f2o = 100, J = 100, a = 1 and different strengths of the disorder, ap = 0; 15; 50. Open symbols stand for the 
numerical data, full symbols, for the analytical estimate (25). 

We turn now to the fidelity. Since in our model there is a large off-diagonal part determined by the constant term f2 
and by random terms ap, we introduce an additional random (small) terms Sp that can be treated as the perturbation. 
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In this way we can analyze the dependence of the fidelity on time t and on the perturbation e'^ = (e^) for different 
values of J and ap. Typical dependencies for the fidelity J-'{t) are reported in Fig. 4. First, it is interesting to discuss 
the data for the dynamical model, when o-p = 0. Unexpectedly, for the integrable case (J = 0) the fidelity is much less 
than for the quantum chaos (J = 100). This fact may have the same origin that was discovered in Ref [36]. However, 
we should stress that our model has no classical limit, therefore, we can speak about integrability or chaos without 
the reference to the classical dynamics. Moreover, in Ref. [36] the effect of a stronger stability for chaotic situation 
occurs mainly on a small time scale. In our case the effect is clearly seen on a large time scale as well. One should not 
also forget that for J = the unperturbed spectrum is highly degenerate, therefore, the observed effect of a stronger 
stability for larger values of J may be related to the break of the degeneracy. 

Another effect is that practically there is no difference for the values ap = 15 and ap = 50. This means that for 
ap > 15 the disorder is quite strong and does not depend on its strength. Indeed, for such a disorder, the difference 
between integrable (J = 0) and chaotic (J = 100) cases is not big. Thus, the results given by the fidelity confirm 
the conclusion drawn from the entropy data: dynamical chaos that is due to a strong interaction J = 100 between 
qubits, and chaos due to a strong disorder cr = 50 are similar (compare fidelity for J = 100, ap = with that for 
J = 0,ap = 50). 




0.5 1 1.5 2 0.5 1 1.5 2 0.5 1 1.5 2 

£t 

FIG. 4. Fidelity for different inter-qubit interactions J and the disorder ap, with f2o = 100. The time dependence is shown 
for 3 different values of the perturbation, e = 5 (squares), e = 10 (circles), and e = 20 (triangles). 

It is also instructive to analyze more carefully the time-dependence of the fidelity for the dynamical case when 
ap = 0, see Fig. 5. One can see a very good scaling !F = T{et) for the integrable case J = 0, in contrast with the 
chaotic case J = 100. Moreover, as is seen from the data in semi-log scale, practically for all values of T, apart from 
very small values, the time dependence is the Gaussian , ~ exp(— Ce^t^). This is in contrast to the chaotic case 
with J = 100 (additional data show that the dependence of T on time is closer to the exponential one rather to the 
Gaussian one). The origin of a good scaling, together with its Gaussian form is not clear, however, one can expect 
that it is related to the form of the strength function, as in the case of the TBRI model. To clear up this problem, 
additional numerical study is needed. 
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1 2 30123456 



FIG. 5. Comparison of the fidelity for the regular, J = 0, and chaotic, J = 100, cases, with f2o = 100. In both cases the 
disorder is absent, Up = 0. Different scales are used for the fidelity and et, in order to reveal the time-dependence on a large 
time scale. As in Fig.4, e = 3; 5; 7 with squares, circles and triangles, respectively. 

V. CONCLUDING REMARKS 

As is mentioned above, there are many studies where the time-dependence of the fidelity is discussed from the 
viewpoint of quantum chaos. The main attention was paid to a long time scale on which an exponential or Gaussian 
decay typically occurs. In the case of the exponential decay the question of general interest is about the characteristic 
parameter of this decay. Specifically, it is widely discussed whether the exponential decrease of the fidelity is governed 
by the Fermi golden rule (in this case the strength function has the Lorentzian form), or by the classical Lyapunov 
exponent of the corresponding classical system. 

On the other hand, in application to the models of quantum computation the problem of a long-time behavior of the 
fidelity seems to be irrelevant. Indeed, for a quantum computation one needs to have a very stable regime where the 
fidelity is close to one. Therefore, in this application one should pay the main attention to the short time scale, where 
the fidelity differs from 1, say, not more than 10~^. This is because the time of the quantum computation can be 
very large, and many pulses of an external magnetic field are needed to implement the quantum protocol. One should 
expect that in this situation the perturbation theory works well (see Ref. [35]). Let us see what the perturbation theory 
gives for our model (32) describing the system within a single pulse. In our case the total Hamiltonian H = Hq + V 
has both diagonal and off-diagonal parts, and let us assume that the perturbation is described by an additional term 
S of the same structure as the fi— term in Eq.(32), with the variance for its matrix elements. Also, we assume 
that the interaction part V may have additional terms, see above. 

In this case the perturbation theory gives, 

jr(f) = 1 - 5%t^ - Re{R)t^ (33) 

where 

5%=Y. El^,^ = e'N, ; R = HE-EH. (34) 
with as the number of non-zero matrix elements of S for a fixed fco • 
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One can see that when the unperturbed Hamiltonian H has a diagonal form (for example, when V ~ 0) the terms 
H and S commute (therefore, R = 0) and we come to the previous expression (13) with — S^. Since in our 
numerical study the perturbation has the same structure as the off-diagonal Q— terms, one can get R = VS — "SV = 
(assuming {T,m,ko) = 0)- Therefore, for our particular type of the perturbation, there is no influence of the term V. 
Our numerical data confirm these findings. The expression (33), however, shows that in a general case the fidelity 
strongly depends on the type of the perturbation. 

In conclusion, we have studied the time dependence of the Shannon entropy and fidelity for the model of a quantum 
computation, in comparison with analytical predictions obtained for the model with two-body random interaction. 
In spite of a big difference between these two models (one is the dynamical one and another is random), we have 
found that in many aspects some properties of the dynamics are quite similar. One of the important results is that 
global properties of the dynamics look the same both for the dynamical model with strong chaos, and for the model 
with a strong disorder. In particular, the entropy and fidelity behave in the same way, manifesting the relaxation of 
the system to a statistical equilibrium. The time scale on which this relaxation occurs can be described by a linear 
increase of the entropy, with a good correspondence with the simple analytical expression. This fact confirms the 
expectation that the TBRI model can serve as the base in understanding the properties of quantum many-body chaos. 
Our numerical data for the fidelity do not confirm the expectation that this quantity may serve as a good indicator 
of the quantum chaos. Specifically, there are many open questions related to the problem of the universal properties 
of the fidelity. Much depends on the type of the perturbation and on the form of an initial packet, therefore, more 
extensive studieds are needed. As for the application of the fidelity to a quantum computation, one can expect that 
real interest is restricted by a small time scale where the standard perturbation theory works very well, and there is 
no infiuence of the quantum chaos. 
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